P(I) 



1. 1 



1.08 



1.06 



1.04 



1.02 



0.2 0.4 0.6 0.8 1 



Figure 1. Periodic Orbits of SAM: Nunes, Casasayas, Tufillaro. 
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Figure 2. Periodic Orbits of SAM: Nunes, Casasayas, Tufillaro 



- (a) 



p= 1.10069 1=1.0 




- (b) 



p = 12/11 1 = 0.8438 




- (c) 



p= 13/12 1 = 0.7343 




- (d) 



p = 14/13 1 = 0.6488 




Figure 3. Periodic Orbits of SAM: Nunes, Casasayas, Tufillaro. 



Periodic orbits of the integrable swinging Atwood's machine 

Ana Nunes 

Departamento de Fisica 

Universidade de Lisboa 

Campo Grande, CI, Piso 4 

1700 Lisboa, Portugal 

Josefina Casasayas 

Departament de Matemdtica Aplicada i Analisi 

Universitat de Barcelona, Gran Via 585 

08071 Barcelona, Spain 

Nicholas Tufillaro* 

Center for Nonlinear Studies and T13, MS-B258 

Los Alamos National Laboratories, Los Alamos, NM 87545 USA 

Submitted to The American Journal of Physics 

Manuscript Number 5708 

Abstract 

We identify all the periodic orbits of the integrable swinging Atwood's ma- 
chine by calculating the rotation number of each orbit on its invariant tori 
in phase space, and also providing explicit formulas for the initial conditions 
needed to generate each orbit. 
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I. INTRODUCTION. 

Integrable Hamiltonian systems typically display an infinity of distinct periodic and 
quasiperiodic orbits. The bounded Kepler problem often studied in classical mechanics is 
not typical in this respect since it only exhibits periodic orbits of a simple type. This is 
because the Kepler problem has a third invariant, the Runge-Lenz vector [1], in addition 
to the energy and angular momentum. This third invariant forces all the bounded motions 
to be periodic. A somewhat more typical example of an integrable two-degree of freedom 
Hamiltonian system is the swinging Atwood's machine when the mass ratio (/x = M/m) of 
the non-swinging to swinging mass is equal to three [2,3]. 

Numerical studies of the integrable swinging Atwood's machine exhibit a plethora of 
distinct periodic and quasiperiodic orbits, and when viewed in configuration space it appears 
to be a difficult task to organize and classify all these different types of orbits. 

The modern view of classical mechanics emphasizes the "geometry of phase space," and 
with this geometric perspective the problem of the classification of orbits in integrable Hamil- 
tionian systems becomes elementary. A basic result of the geometric approach to classical 
mechanics states that the topology of the invariant manifold M of integrable Hamiltonian 
systems must be an n-dimensional torus [4]. In particular, for a two-degree of freedom 
Hamiltonian system with just two invariants, each orbit in the four-dimensional phase space 
must be confined to a two-dimensional torus, a hollow donut. Furthermore, all orbits which 
arise from flows (vector flelds) on a two-dimensional torus must be either a periodic or 
quasiperiodic winding of the torus. In the latter case the quasiperiodic orbit will densely 
cover the whole torus. 

The "natural" way then to classify the orbits is to calculate rotation number about the 
torus: the ratio of winding about each generator of the torus. For quasiperiodic orbits this 
rotation number is irrational, and it is a single number that uniquely identifles the orbit. 
In the case of a periodic orbit, this number is rational and it uniquely identifles a family of 
periodic orbits all of the same type or shape when viewed in conflguration space. Thus the 



rotation number, which uniquely indexes each torus (and hence each type of orbit) in the 
phase space of an integrable Hamiltonian system, provides the most natural solution to the 
problem of identifying the orbits of a two-degree of freedom integrable Hamiltonian system. 
In this paper we use some classical results of classical mechanics (Hamilton- Jacobi theory) 
along with the modern perspective of classical mechanics (the geometry of phase space) in 
order to identify and classify all the periodic orbits of the integrable swinging Atwood's 
machine. It is hoped that this elementary but physically realistic example might find use in 
an advanced classical mechanics course which attempts to introduce the flavor of a modern 
approach to classical mechanics. 

II. THE ORBIT EQUATION. 

Taking polar coordinates in the vertical plane where the weight of swinging mass m 
moves, the total energy of the system is 



E = T + V 

-r^ + -mr^e^ + gr{M - m cos 9). 



M + m .^ 1 2/12 



The conjugate momenta pr and pe are 



Pr = {M + m)r, (2a) 

Pe = mr^e. (2b) 

So the Hamiltonian of the problem is given by 

H{r,e,pr,p8) = 

- ( ^ — + ^]+gr(M -m cos 6). (3) 

As shown in [3], it is convenient to change to another set of canonical coordinates where the 
equations of motion become simpler, especially if we want to do Hamilton- Jacobi theory. 
Let the coordinate transformation to (^,77) be deflned by 



( = ^r{l + sin(^/2)), 
rj = ±^r{l - sin(^/2)). 



(4a) 
(4b) 



The inverse change of coordinates is 



e = 2arctan [{(^ -ri^)/2(ri] 



(5a) 
(5b) 



and the generating function F2{r,6,p^,pr^) which determines the extension of this change 
of position coordinates to a canonical transformation in phase space will be 



F2{r,0,p^,Pr,] 



^r(l + sin(^/2)) p^ + ^r(l - sin(^/2)) p^. 



(6) 



Hence, 



Pr 



dF2 
dr 



^1 + sin(^/2) ^1 - sin(^/2) 



2v^ 



-Vi + 



2V^ 



Pv 



Pe 



e+v 



dFo 



(7a) 



36 



r cos 



(9/2) 



--Pi 



r cos 



{612) 



4^1 + sin(^/2) 4^1 - sin(^/2) 

VPi - CPv 



Pv 



(7b) 



Since we know (see [5]) that the system is integrable only when M/m = 3, we set M = 3 
and m = 1 throughout this paper. For these values of the parameters, eq. (3) becomes, after 
introducing eqs. (5) and eqs. (7), 



HU,V,Pi,Pv] 



e+v' 



lipl+pD + Ht + v' 



(8) 



the expression of the Hamiltonian of the system in the new set of coordinates. 

We recall that the idea behind Hamilton- Jacobi theory is to find a generating func- 
tion S'(^,77,/3i,/32, t) for a canonical change from {( , f] ^ Pi ^ Pv) ^° ^ ^^^ ^^^ °^ variables 
(ai, q;2,/3i,/32) which are constant along the motion, i.e., such that the Hamiltonian ex- 
pressed in terms of the new variables is identically zero. The function S with this property 
satisfies the first order partial differential equation 

„/, dS dS\ dS ^ 

and the constants /3i,/32 will be the two independent constants which will come out of the 
integration of eq. (9). Therefore, by construction, solving eq. (9) is equivalent to solving the 
equations of motion. 

As shown in [6], eq. (9) for Hamiltonian eq. (8) separates and yields a solution in terms of 
quadratures, i.e., the problem reduces to that of computing certain integrals. More precisely, 
let us consider eq. (9), where H is the Hamiltonian given by eq. (8), and try a solution of 
the form 

S{(,V,t) = S^{0 + S2{v) + Ss{t). (10) 

Then, eq. (9) becomes 



1 



e+v' 



Kf^t>^'«*^'^' 



Clearly, each member must be equal to a constant, /3i, and moreover, since by definition 

"" (12a) 

(12b) 

this constant must be the total energy of the system, i.e., 

/3i = E, (13) 
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Pi = 


dC 




dS 


Pv = 


drj' 



ana so 



Ss{t) = -Et. (14) 



Equation (9) then reduces to another separable equation in ^ and 77, 



2 



1 / db\ 1 , r. f4 77i(i-2 



8 V ^^ 



+ ^gc - Ee 



2 



1 / 0,02 \ r. 4 I 771 2 



2577^ + Erj'. (15) 



8 \dri J 

Again, both members must be equal to the second integration constant /32, which we shall 
name 1/8. Therefore, denoting by s{x,E,I), 



s[x 



E, I) = r J I - 16gy^ + 8EyHy, (16) 

Jk 



where k is some constant to be fixed later, the solution S[(,r], E, /, t) is, up to two indepen- 
dent additive constants, 

5(^, 77, E, /, t) = s{C, E, I) + s{rj, E, -I) - Et. (17) 

The equations of the canonical transformation whose generating function is S which will 
give the complete solution of the problem are 

dS ds ,^ ^ ^. ds , ^ ^. , . 

a2 = -^^{C,E,I) + -^^{rj,E,-I), (18b) 

P^ = ^ = |(^,^,^), (18c) 

P. = |^ = |^('7,^,-/)- (18d) 

Note that the two independent additive constants are now included in eqs. (18a) and (18b). 
Equation (18b) is the orbit equation, which gives a relation between the position coordinates 
( and T] in terms of the constants E, /, and 0:2. Adding the information contained in eq. 
(18a), that contains time explicitly, we shall obtain the trajectories equations. It is therefore 
clear that the other constant ai plays the role of the initial time. The last two equations of 
(18) complete the description of the motion in phase space. 



Let us now focus our attention on the orbit equation, to obtain an explicit expression in 
terms of known functions. First, we note that, as pointed out in [7], we may always choose 
appropriate time and length units to set E = g = 1. From now on, we shall assume that 
we are working in these units. In particular, the parameter / will always take values in the 
interval [—1,1]. Then, the orbit equation becomes 

1 /^ dy 1 /'' dy 



"2 



/■« dy i r^ 

Jk, JI - 162/4 +1P ~ 2 A, J^ 



2 h, ^I - 162/4 + 82/2 2 A, ^-I - 162/4 + Sy^ 
= ^S(^,/) + ^0('7,/)- (19) 

Fixing the constants h\ and h2 at appropriate values and assuming for now that / > 0, we 
obtain 

^{U) = Y^{U\ (20a) 

0(77,/) = ^2(77,/), (20b) 



where 



ri(x, /) = r , ^^ (21a) 

r2(x, /) = r , "^^ (21b) 

h ^(c2-2/2)(2/2-d2) 



anc 
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1 + VI + / 


.2 


4 
1 - Vl - / 



6^ 



2 -i + ^TTT 



4 

On the other hand, we also have, for / < 0, 



4 



^2 i + vr^7 



S(^,/) = -r2(^,-/), (23a) 

Q{rj,I) = -Y,{rj,-I). (23b) 

Therefore, the orbit equation may always be expressed in terms of the functions Yi[x,I) 
and 1^2(33,/), / > 0. Now, both Yi and Y2 are elliptic integrals of the first kind (see, for 
instance, [8]). More precisely, we have, denoting as usual by F[ip\m) the elliptic integral 



F((p\m)= rCi-msin^ (l))-^l^d(l) 
Jo 

/•z=sin if 

= / [{l-y'){l-my')]-'/'dy, 

Jo 



(24) 



wi 



ith 



anc 



Yi{x,I) 
Y2{x,I) 



ViTl 



^(t/5i|mi^ 



1 + ^r^T 



F{ip2\m2] 



mr{I) 
MI) 



1 + ViTl 

2VTn ' 
2^r^7 



i + VT^' 



sin(^i(x, /) 



8x^^/lTl 



\ 4x2(1 + ^rTT) + /' 



1 - 4x2 + ^r^7 



\ 2^r^7 



s'mip2{x,I) 

Hence, for / > the orbit equation may be written 

:F(c^i(^,/)|mi(/)) 



1 



"2 



V2T2/ 



+ - 



V2 



,F[^2[vJ)\m2[I)l 



(25a) 
(25b) 



(26a) 
(26b) 



(27a) 
(27b) 



(28) 



"^|l + v^^T" 

where (^i(^,/), TOi(/), ip2{r],I), and 7722(1) are given by eqs. (26) and (27). For / < 0, by 
eq. (23), the orbit equation is obtained from eq. (28) by interchanging ( and rj. 

Using eq. (28) and a symbohc math program, and going back to the original coordinates 
r, 6, it is easy to obtain plots of the orbits of the system in the plane. 



III. LOW-ORDER PERIODIC ORBITS (TORUS KNOTS) 

The system we are dealing with is non-degenerate and so we expect to find periodic and 
quasi-periodic orbits, according to the value of / that we pick in a certain energy level. Since 
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we are studying an integrable Hamiltonian system, all these orbits, both the periodic and 
the quasi-periodic will lie on tori in the phase space, which form the closure of the quasi- 
periodic orbits in one case, and which will be foliated by the periodic orbits in the other 
case. This determines the knot type of the periodic orbits: all of them must be torus knots. 
A torus knot is characterized by two integers, n,m, that measure the number of times the 
curve winds around each generator of the torus. The rational rotation number p associated 
to each of the tori foliated by periodic orbits will measure the ratio n/m. Therefore, if we 
want to know which torus knots are realized as orbits of our system, or if we want to know 
how to pick the "simplest" torus knots, i.e., those for which both n and m are small integers, 
we need information on the dependence of the rotation number on the parameter /. The 
range of the function /)(/) will tell us which are the low-order resonances and hence which 
are the low-order torus knots, and the form of the function itself will tell us where to look for 
these knots, i.e., which approximate value for / we should take. The appropriate variables 
to study this question are action-angle variables. 

Recall that, in general, the action variable Ji is given by §pidqi, where the integral is 
taken along a period of the projection of the orbit on the qi,pi plane. The angles are the 
canonical conjugate of the actions, and the Hamiltonian in this set of variables depends 
only on the Ji. The angles change with time with constant frequencies, given by the partial 
derivatives of H with respect to the Ji. Hence, in our case, the rotation number associated 
to the motion on the two-dimensional tori that foliate the energy level E = 1 will be given 

by 

where the derivatives must be evaluated along H[Ji, J2) = 1. Now, from equation (6.3) and 
(6.4) we have 

= j ^I + 8Ee - IQC'dC = j{E, I) (30a) 



J2{E, I)= j Pr,dr] = j -K-dv 

= i ^J-I + 8^772 - IGrj^drj = j{E, -I), (30b) 

where we have recovered the dependence on E since it is essential for the computation of 
eq. (29). From the identity E{Ji{E, I), J2{E, I)) = E we get 



dE/dJ, ^ '-ifjlj) ^ i(l,-/) 
dEldJ, ^{1,1) §2(1,/) 



Pin = ^^7^ = -#7T^ = ^iTTf^- (31) 



But 



i(^.^)4/^TTlF^ (-' 



and so, for / > 0, 



/ 



dy 



^(c=-y=)(y=-d=) 



Pil) = r ay (33) 



where a, h, c and c? are given by eq. (22), while as in eq. (23), 

p{-I) = l/p{I) (34) 

holds. Hence, it is enough to study the function p on the interval [0, 1]. This last identity 
arises from the time reversibility of the system since the first integral is invariant under a 
change in the velocity sign. 

As we have already seen, the indefinite integrals correspond to elliptic functions of the 
first kind. Now, the integrals in eq. (33) are the corresponding complete elliptic integrals of 
the first kind K[m) (see, for instance, [8]). In terms of these functions, eq. (33) becomes 



Pin 



2^1 + 7 K[m,{I)] ^g^^ 



\l + ,/l^K[rm{I)y 
where mi[I) and 7722(1) are given by eq. (26). This expression and the range |o([0,l]) can 
be easily computed using for K[m) the polynomial approximation 

K[l — m) = (flo + a-iTn + a2m^) + 

{bo + bim + b2m'^)log{l/m) + e{m), (36) 
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where |e(m)| < 3.10"^ and ao = 1.38629, a^ = 0.11197, as = 0.07252, bo = 0.5, b^ = 0.12134 
and 62 = 0.02887. 

Figure 1 is a plot of p[I) verses /. The graph shows a monotonic function with p[0) = 1 
and p{l) = 7r/(2^/^iC[(l + y/2)/2^^^]) % 1.10069. Thus, no periodic orbit can exist in the 
integrable swinging Atwood's machine with rotation number greater than about 1.10069. 
The "simplest" periodic orbits will therefore be of fairly high period and will have rotation 
numbers like 11/10, 12/11, 13/12 and so on. 

In order to plot some of these periodic orbits it is useful to express the initial conditions 
as a function of /. Then to locate a periodic orbit with a given rational rotation number, 
all we need to do is calculate the inverse function p~^[n/m) = I and then use this value of 
/ to find the initial conditions needed to generated an orbit with rotation number p = n/m. 

To find a formula for the initial conditions we note that at the surface of section ^ = 0, 
the energy E{r,6,r,6) and first integral I[[r,6,r,6) are with E = g = 1 [9]: 

2r^ + -rH^ + 2r = 1, (37) 

and, 

/ = 16rV^. (38) 

Eliminating 6 between these two invariants we find 

1 



r[r. 



I) = ±^> 



--r] ± 



1 -' 
r 



(39) 



V2^|V2 J \\2 J (16r)2- 

All real solutions of eq. (39) with / G [0, 1] represent cross-sections of the invariant tori at 
^ = 0. A few of these cross-sections are shown in Figure 2. The "first" torus occurs at 
/ = and corresponds to the family of "teardrop-heart" ejection-collision orbits previously 
studied [6]. The "last" torus occurs at / = 1 and corresponds with the "loop" periodic orbit 
shown in Figure 3(a). This last torus is degenerate, and it is the only periodic orbit in the 
integrable swinging Atwood's machine which may have a non-rational rotation number. All 
other orbits on the tori between < / < 1 are either a family of periodic orbits if p[I) is 

11 



rational, or a single quasiperiodic orbit otherwise. Equations (39) and (38) with r = 0.25 
and ^ = provide explicit initial conditions for an orbit with any value of /. The "last 
torus" (loop orbit) has the initial conditions (ro, ^o,ro, ^o) ~ (0.25,0,0.353553,2.82843). 

Figure 3 illustrates a few of the simpler periodic orbits that occur in the integrable 
swinging Atwood's machine. 

IV. LOOP ORBIT 

We conclude by showing that it is possible to write down an explicit exact solution for 
the loop orbit [E = g = 1). In the (^,77) coordinates, initial conditions for the loop orbit 
are C = V — V^ with ^ = v2/2 and rj = 0. And it is easy to check that 77 = 0. So, the 
solution of the loop orbit in the (^,77) coordinates is particularly simple, namely: ( = ([t) 
and T] = 1/2 for all t. Now, using eq. (5), we find that on transforming back to the original 
(r, 6) coordinates this orbit gets mapped to: 

r_{e) = \ sec(f ) [sec(f ) - tan(f )] for < ^ < 27r, (40a) 

r+(^) = |sec(f)[sec(f) + tan(f)] for - ^0 < ^ < ^0 (40b) 

where the range of the solution on the second branch of the orbit (eq. 40b) can be determined 
by setting the velocities {r.,0) equal to zero in the energy equation (1), 

r[3 - cos(^)] = 1 (41) 

and the second invariant I{r,6,r,6) [9], 

16r^in (^^^ cos^ (^^^ =1, (42) 

which is equal to one for the loop orbit. Solving for 6 we find: 

(9o = arccos(-5 + 4^2) ^ 0.854 radians. (43) 
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Figure Captions 
Figure 1. Rotation number p[I) as a function of the invariant /. 
Figure 2. Cross sections (at ^ = 0) of invariant tori for a few values of /. / = 

corresponds to the "teardrop-heart" orbits while 1=1 corresponds to the "loop" orbit, a 

degenerate torus at r = 0.25. 

Figure 3. A few periodic orbits and their rotation numbers with initial conditions 

(ro,^o,ro,^o): (a) (0.25, 0, 0.353353, 2.82843), (b) (0.25, 0, 0.438269, 1.92533), (c) (0.25, 0, 

0.458094, 1.603), (d) (0.25, 0, 0.469171, 1.38283). 
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